The Magnetic Ordering of the 3d Wigner Crystal 
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Using Path Integral Monte Carlo, we have calculated exchange frequencies as electrons undergo 
ring exchanges of 2, 3 and 4 electrons in a "clean" 3d Wigner crystal (bcc lattice) as a function of 
density. We find pair exchange dominates and estimate the critical temperature for the transition 
to antiferromagnetic ordering to be roughly 1 x 10~*Ry at melting. In contrast to the situation in 
2d, the 3d Wigner crystal is different from the solid bcc ^He in that the pair exchange dominates 
because of the softer interparticle potential. We discuss implications for the magnetic phase diagram 
of the electron gas. 



I. INTRODUCTION 

The uniform system of electrons is one of the basic 
models of condensed matter physics. Wignei»l pointed 
out that at low density, the potential energy dominates 
and the system will form what is now called a Wigner 
crystal (3dWC). There have been attempts to make lab- 
oratory examples of the low density homogeneous 3d elec- 
tron gas with specially designed band-engineered AlGaAs 
heterostructures^. In this paper, we report on calcu- 
lations of the spin Hamiltonian in the low density 3d 
Wigner crystal. 

At r = 0, the properties of the electron gas are deter- 
mined by a single dimensionless parameter = a/ao = 
{m* /me) a /at, where a = (3/47rp)^/'^, ab is the Bohr ra- 
dius, TO* is the effective mass, and e the dielectric con- 
stant. In this paper we will use effective Rydbergs for 
energies Ry* = (TO*/mee^)Ry and a for units of length. 
In these units the Hamiltonian is: 
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Though a variety of methods have been applied to cal- 
culate the properties of the low density electron system, 
the most successful have been direct simulation methods: 
namely Quantum Monte Carlo. Ceperley and Alder^ us- 
ing Diffusion Monte Carlo (DMC) determined that melt- 
ing at zero temperature occurs at ~ 100 ± 20 for spin 
1/2 fermions, and at ~ 160 for bosons. (This is also 
the zero temperature melting density for distinguishable 
particles). Later estimates^ found melting at higher den- 
sities = 65 ± 10), but the variational trial functions 
in the liquid phase were not sufficiently accurate. Very 
recently, Drummond et al& confirmed the estimate of 
rs = 106 ± 1 using DMC with a variety of better func- 



tions and the more accurate results of Zong et al^ in 
the fluid phase. The bcc crystal structure has the low- 
est energy throughout the stability region of the crystal. 
See, for example, the QMC calculations of Harris et alA 
who found bcc the lowest energy structure in the range 
> 60. 

Once the melting density is established, it is of inter- 
est to determine the low temperature spin order. Har- 
ris et alA reported finding the ferromagnetic bcc phase 
as stable, however, other aspects of those calculations 
have not been reproduced. Drummond et alA attempted 
to determine directly the energy difference between fer- 
romagnetic and antiferromagnetic orderings using DMC 
but found the difference zero within their error bars (on 
the order of 2 x 10"'' Ry at = 100). This is con- 
sistent with the results reported below. One needs to 
use a method sensitive to the small magnetic energies, 
which are typically many orders of magnitude smaller 
than the plasmon energies that determine the accuracy 
of the DMC energies. 

Jones and Ceperleyi studied the quantum melting 
curve for distinguishable particles. At densities for > 
100 the melting is classical, and occurs for temperatures^ 
kBTmeit = 2/{Tcrs)^y where « 173. We only use 
this melting to determine the region where the crystal is 
stable, and thus the region where magnetic ordering is 
relevant. Fig. ^ summarizes the 3deg phase diagram. 

For spin 1/2 particles, the magnetic ordering is not 
fixed by the spatial ordering. The earliest quantitative 
calculation was using a Slater determinant of localized 
Wannier functions (i. e. Gaussians) by Carr.^ He found 
an antiferromagnetic phase in the Wigner crystal at in- 
termediate density and ferromagnetic at lower density. 
We will compare with this calculation later in the paper. 
Edwards and HillelAS, using a Hartree-Fock method with 
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FIG. 1: Phase diagram of the 3d electron gas showing the 
region of stabiHty of the crystal--, the polarization transition 
from QMC calculations- and the antiferromagnetic transition 
(this work), (color online) 
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FIG. 2; The pair correlation function for the 3dWC (solid 
line) [ts = 100) and for bcc ^He (dashed red line)also at 
melting (24.23cm'^/moie. The functions are nearly identical, 
though the electrons can get significantly closer together than 
helium atoms can. (color online) 



a flexible, delocalized basis and a variety of spin orderings 
found antiferromagnetic ordering at intermediate density 
10 < < 40 and ferromagnetic at lower densities. But, 
as mentioned above, the crystal phase is only stable for 
r, > 106. 

Herringii reviewed the situation as understood in 1965 
in some detail, including the contribution of ring ex- 
changes of electrons. Thoulessi^ introduced the current 
theory of magnetism in quantum crystals. According 
to this theory, in the absence of point defects, at low 
temperatures the electrons will almost always be near a 
lattice site. If the system is constrained to stay in the 
neighborhood of the two perfect lattice positions Z and 
PZ where P is a permutation of particle labels, the ex- 
change frequency equals the splitting between the anti- 
symmetric spatial state and the symmetric spatial state: 
2Jp = Ea — Es > 0. The spin Hamiltonian comes about 
from making the total wavefunction antisymmetric: 

'^spin ^ ^ ( 1) Jp^spin (2) 

P 

where the sum is over all cyclic (ring) exchanges de- 
scribed by a cyclic permutation P, and Pspin is the corre- 
sponding spin exchange operator. (Although more com- 
plex products of several ring exchanges are possible, in 
cases considered, they are negligible.) The sign, (— 1)"^, 
implies that an exchange of even number of electrons is 
antiferromagnetic and an odd number of electrons is fer- 
romagnetic. Ring exchange models have been used to 
describe correlated electron systems, such as high tem- 
perature superconductors^"*, quantum Hall systemsl4 as 
well as electrons^^ and helium atoms^^ confined in planes. 



One might expect that pair electron exchanges would 
dominate over higher-body exchanges. Since the bcc lat- 
tice is bipartite, a simple Neel antiferromagnetic state 
would seem to be favored. Rather surprisingly, it has 
been foundi^i that in 3d solid '^He, which also forms a 
bcc lattice, exchanges of 2, 3 and 4 particles have roughly 
the same order of magnitude and must all be taken into 
account to understand the magnetic ordering. This is 
known as the multiple spin exchange model(MSE). The 
resulting spin order is more complex since the order is 
frustrated by the competing exchanges. We wish to de- 
termine whether such a model is relevant for the 3dWC. 

Figure 121 shows the pair correlation functions for solid 
■^He and the electron gas near the crystallization density. 
Because the (?(r)'s are so similar, one might expect their 
exchanges frequencies would be similar and hence have 
the same magnetic ordering. We also note that the Lin- 
demann's ratio, the mean squared displacement in units 
of the nearest neighbor spacing, for bulk helium and the 
Wigner crystal are also similar near melting (0.32 and 
0.30 respectively). However, note the g{r)'s are very dif- 
ferent at small r because the potentials are so different; 
the helium-helium interaction is much more repulsive at 
short distances. 

In this paper, we determine the magnetic interaction 
in the Wigner crystal, based on Thouless'^^ theory of 
exchange. Path Integral Monte Carlo (PIMC) as sug- 
gested by Thoulessi^ and Roge»i^ has proved to be a 
reliable way to calculate these parameters directly from 
the Coulomb interaction. The theory and computational 
method have been tested thoroughly on the magnetic 
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properties of bulk helium obtaining agreement with mea- 
sured propertieaiSi. Wei^ have also used this method to 
calculate exchange frequencies in the 2dWC. Here we re- 
port results for the 3dWC. 



II. COMPUTATIONAL DETAILS 

The Path Integral method for calculating exchange fre- 
quencies is based on the ratio: 



/p(/3) 
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(3) 



where Z represents the many-body configuration of elec- 
trons sitting on the bcc lattice sites and PZ is a permu- 
tation of those sites. Then under general assumptions, 
the exchange frequencies are given by: 



/p(/?) =.tanh(Jp(/3-/3o)) 



(4) 



where (3q is the amount of imaginary time to initiate the 
exchange. The ratio fp is determined by a specialized 
Path Integral Monte Carlo method and Eq. 2|is inverted 
to determine Jp. 

We do simulations with two types of paths: a) paths 
beginning and ending at the perfect lattice positions and 
b) paths beginning at Z and ending at PZ. The imagi- 
nary time density matrices e~^^ are expanded out into 
a path integral connecting the end-points of the paths. 
Using the polymer "isomorphism"—, fp{z) is related to 
the free energy need to induce a specific cross-linking 
P into a crystal of ring polymers. Using the Bennett 
methodSSiSi, we directly determine the ratio fp by ex- 
amining histograms of the change in action in mapping 
paths of one type into paths of the other type. With this 
method we can determine very small frequencies to an 
accuracy of a few percent, irrespective of the magnitude 
of Jp. 

Ewald sums are used to represent the Coulomb inter- 
action in periodic boundary conditions, taken as a cube. 
The potential is split into long-range and short-range 
terms with the usual Gaussian breakup^"^. Then the ex- 
act two particle action for the short ranged potential (the 
complementary error function) is determined using ma- 
trix squaring;™ The long range action is taken in the 
primitive approximation; this is appropriate since it is 
smooth. Most calculations were done with 54 electrons 
in the simulation cube with a few tests of 128 electrons 
giving agreement. Note that at low density, long wave- 
length plasmons are strongly suppressed by the poten- 
tial energy, making the exchange more localized spatially 
than is the case with solid helium. 

One numerical approximation is the imaginary-time 
step, or equivalently, the number of points on the path. 
We did several calculations for each value of and each 
type of exchange to establish that the results are in the 
zero time step limit within error bars. A second approxi- 
mation concerns the value of P in Eq. 13 Because the ex- 
changes are instantons (i e. confined in imaginary time). 



the results converge quickly in (3 as long as /3o < (3. We 
typically choose P < 3/3o and observe very weak depen- 
dance on /3. Typically, for converged results, this implies 
one to two hundred steps in the exchange. Further details 
of the method have been discussed in earlier papers^ 

In order to keep the system from melting for densi- 
ties higher than the bosonic melting'! 100 < rg < 160, 
the electron paths are restricted to stay in the positive 
region of a trial wavefunction: the fixed-node boundary 
conditions. We separate the electrons into those exchang- 
ing (say "p" of them) and the spectators { i. e. N — p 
electrons). A Slater determinant of only the spectator 
electrons is constructed: 



4'{R(t)) = det[exp(-c(r,(t) - 



(5) 



where Zj is the set of spectator lattice sites and it) the 
imaginary time path of the spectator electrons. The value 
of the parameter c was optimized by variational Monte 
Carlc^^: c — 0.2ry^. The more recent values obtained 

p. 1/2 

by minimizing the fixed- node energy c = 0.1 Irs were 
not available when we did the calculations. A Slater de- 
terminant of Gaussian orbitals is an accurate represen- 
tation for the nodes of the many-electron ground state 
wavefunction. More complicated forms, such as linear 
combinations of Gaussians, do not result in a significantly 
better trial functions^ for the 3dWC. 

Only paths which keep a positive determinant through- 
out the path are kept il}(t) > for < t < /3: this is the 
fixed-node method. We find that these boundary condi- 
tions are sufhcient to keep the system from melting. We 
do see a suppression of the exchange frequency caused 
by the determinantal boundary conditions at Vs = 150 
of about 10%. In principle, the spin ordering should be 
determined self consistently. For example, it would be 
better to apply antiferromagnetic boundary conditions, 
however, we have not tested that approach. For this rea- 
son, we may have corrections to the exchange frequencies 
on the order of 10%, particularly near melting. 



III. EXCHANGE FREQUENCIES 

We have calculated 2, 3 and 4 particle exchanges for 
six different densities in the range 100 < < 150. We 
determined two different pair exchanges: first neighbor, 
Jnm and second neighbor, J„„„. We do find a significant 
2 body next-nearest-neighbor exchange. The most com- 
pact three-body exchange, Jt, has two first-neighbor and 
one second-neighbor electrons. There are two types of 
four electron exchanges involving only first neighborsii: 
the planar exchange, Jp, and the folded exchange, J/. 
Table I shows the exchange frequencies of our calcula- 
tions versus density. 

Figure|31shows the exchange frequencies versus density. 
We find that the exchange energies are very much smaller 
than the plasmon energies. The kinetic energy of the 
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Wigner crystal can be expanded in a power series^A: 

T=-^^^^Q.&mr-^'^ -Q.bbZr-\ (6) 

In the density range of consideration, the kinetic energy 
varies from 0.6 mRy to 0.3 mRy: five orders of magnitude 
greater than the exchange frequencies. Roughly speak- 
ing, the electrons vibrate around their lattice sites 10'^ 
times before exchanging with a nearby electron. This is 
comparable with the situation in bcc ^He and justifies 
that we can reduce the original Hamiltonian involving 
charges, Eq. to the spin Hamiltonian, Eq. |2| 

We see in Fig. O that the exchange frequencies drop 

1 /2 

off exponentially: roughly as exp(— Sr., ). This follows 
from assuming a single most probable path for the ring 
exchangsi^ is independent of density. (Note that we do 
not make this assumption in the PIMC simulation). It 
is likely that WKB calculation of S will give reasonable 
estimates of the exchange frequencies as they do for the 
2dWGi2i^. 

We find that the rate for pair exchanges is much larger 
than the other exchanges, making a very stable antifer- 
romagnetic ground state. The planar 4-body exchange, 
also antiferromagnetic, is slightly larger than the ferro- 
magnetic, 3-body exchange. Only ratios of the exchange 
rates can enter into determining the stability of a given 
magnetic state. Figure 0] shows the density dependance 
of the ratios Jf/Jn and Jt/J„. We see a decrease in 
these ratios as density decreases, thereby further stabi- 
lizing the antiferromagnetic phase. Next in importance 
is the 2-body next nearest neighbor exchange, followed 
by the folded 4-body exchange. 



IV. MAGNETIC PROPERTIES 

Having determined the exchange frequencies, we can 
now discuss the magnetic properties, based on the spin 
Hamiltonian in Eq. |21 In principle, one has a formidable 
many-body problem. However, we can make use of the 
extensive results available from studies of solid ^He, also a 
bcc lattice. In particular, see the review of Roger et alm^ 
(RHD) and references therein. For spin 1/2 particles, 
one can write the spin Hamiltonian in terms of Pauli 
spin matrices. Two and 3 particle exchanges map into a 
Heisenberg model: 

. (1) . (2) 

• C^J + • • ■ (7) 

i,j id 

where the sums are over first, second and third neighbor 
pairs respectively. (Note on notation: Jx > refers to a 
ring exchange frequency for the cycle x, while j„ refers to 
the coupling constant in Eq. between two spins a dis- 
tance X apart; all frequencies have an opposite sign from 
the notation of RHD.) These spin couplings are given in 
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FIG. 3: (color online) Exchange frequencies (in Ry log scale) 
versus r]''^ . The solid (black) line shows the 2-body nearest 
neighbor exchange. The dashed (red) line is from the calcu- 
lations of Carj-. The other lines are the triple (blue short 
dashes) and planar 4-electron (cyan, long dashes) exchange 
frequencies. A WKB calculation would be a straight line. 

terms of the ring exchanges by: 

Jl = Jnn+H~2Jt + Jp + Jf) (8) 
j2 ~ 2( — 2J( + Jf) + Jp + Jnnn 
J3 - Jp/2. 

The pair couplings, shown in Table I, can either be pos- 
itive (antiferromagnetic) or negative (ferromagnetic) de- 
pending on the relative importance of even and odd ring 
exchanges. 

The four spin exchanges lead to additional four spin 
terms in the effective Hamiltonian in Eq. \7\ 

Gi-jki ~ i<7i-crj){ak-<Ji) + {ai-ai){aj-ak) — {(yi-crk){(7j-cri) 

(9) 

The summation (x — p oi x — /) is over distinct 
labels describing the planar or folded four-particle ex- 
changes. According to RHD, one can neglect these ad- 
ditional terms in estimating properties at high tempera- 
tures. Since they are fourth order in the order parameter 
field, they can only contribute at lower temperature. 

Given the above spin Hamiltonian, several things can 
be easily computed. At high temperature, the Curie- 
Weiss constant measures the leading term in a 1/T ex- 
pansion to the susceptibility (x^^ = C{T-Q+B/T+. . .). 
We find that: 

e = Unn + 3J„„„ - 36Jt + 18Kf + 18Kp. (11) 
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TABLE I: Results of the PIMC calculations at various values of r^. Units are 10"^ Ry/electron. Quantities in () are the 
estimated relative statistical error. The second column is the estimated transition temperature (also in nRy) based on Eq. 1131 
The spin couplings j„ are defined in Eq. |H| 
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FIG. 4: (color online) Ratio of exchange frequencies versus 
Vs. Shown are Jt/ Jnn (blue, short dashes) and Jp/Jnn (cyan, 
long dashes). 

In regions above the ordering temperature, this deter- 
mines the deviation with respect to uncoupled spins. It 
is positive (antiferromagnetic) , for all densities. 

From the dominance of Jnnn over the other ring ex- 
changes, we expect ordering into an antiferromagnetic 
phase via a second order transition at a sufficiently low 
temperature. The complications resulting from a frus- 
trated spin Hamiltonian are absent, so that we can be 
considerably more confident in our predictions than is 
the case with a frustrated spin-model such as helium or 
the 2dWC. First, let us determine the transition temper- 
ature within mean field theory. RHDi"^ (Fig. 16) show 
the mean field phase diagram as a function of (ji, j2i Ja) 
As mentioned above, the four spin terms will not change 
this diagram. Based on the values in Table I, we find that 
the 3dWC is very much in the antiferromagnetic region. 
The density dependance shown in Fig. 0] has a tendency 
to enhance even more the antiferromagnetism at larger 
values of Tj. 

Within mean field theory, the antiferromagnetic phase 



has a critical transition temperature: 

keTc = 4ji - 3j2 - 6j3. (12) 

We can include fluctuations into this estimate by using 
results of series expansions. Oitmaa and ZhengS^ have 
calculated properties of the ji , j2 model on a bcc lattice 
using Fade approximates to high temperature expansions 
(up to tenth order in the inverse temperature.) They find 
the antiferromagnetic phase is stable at low temperatures 
for j2 < 0.705ji. They also find that the transition tem- 
perature is ksTc « 2.76ji — 2.61^2. The effect of fluctu- 
ations is to lower Tc by about 30%. Since the effect of j2 
and js is to couple the same sublattices, a reasonable way 
to extend these results to 7^ is to assume that the 
effect of j3 is determined by the number of spins coupled. 
Since there are twice as many third neighbors as second 
neighbors, we assume: 

fcsT, ~ 2.76J1 - 2.61(j2 + 2j3). (13) 

The lower line on Fig. ^and the second column of Table 
I show this estimated transition temperature. 

For = 120, where we have done the most cal- 
culations, we find an estimated transition temperature 
of Tc = 2.0 nRy. We can parameterize the density 
dependance of the transition temperature by assuming 
the WKB form for the dependance of the exchanges 

1 /2 

on the density. We obtain ksTc « 3.7exp(— 1.92rs ) 
Ry. (Note: because we calculated J„„„ and J/ only at 
Ts — 120, we estimated their values at other densities by 
scaling. The importance of these exchanges is negligible.) 

We now compare the present calculation with DMC 
calculations aimed at determining the spin ordering. This 
is done by performing a ground state total energy calcu- 
lation (within the fixed node method) for a fully polar- 
ized Wigner crystal and for an antiferromagnetic crystal. 
Such an attempt was made both for the 2dWC?§. and for 
the 3dWG^ but without finding a significant energy dif- 
ference. Given the exchange frequencies we can estimate 
this energy difference. Assuming classical spins, the anti- 
ferromagnetic and ferromagnetic energies are easily found 
to be: 

Eaf/N = -2ji-f 1.5j2 + 3j3 (14) 
Ef/N = 2ji-f 1.5j2 + 3j3. (15) 
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Oitmaa and Zheng^^ determine corrections to the clas- 
sical antiferromagnetic ground state energy of the (^1,^2) 
model using an expansion technique and find: 

EafIN = -2.3ji + 1.32J2. (16) 

This is only roughly 10% lower than the classical value 
for the 3dWC parameters. Note that the classical fer- 
romagnetic energy needs no correction. Using these es- 
timates and Eq. |S1 within a few percent, the energy 
difference to spin polarize the crystal is: 

« {Ep - Eaf)/N = 4.3ii + 0.18(j2 + 2^3). (17) 

Near the melting density of = 100, A£; « 2.4 x 10"® 
Ry which is a factor of 10 smaller than the reported^ sta- 
tistical errors within DMC. Even if one could reduce the 
statistical error, the usual spin wavefunction in the anti- 
ferromagnetic phase is relatively crude, and the nodal 
surfaces are not necessarily optimized for the regions 
where exchanges occurs. Calculation of the ring exchange 
frequencies are a much more direct and efficient way of 
determining the magnetic ordering, and yield more phys- 
ical insight into the microscopic mechanisms giving rise 
to the magnetism. 

Concerning previous calculations on the magnetic 
properties of the bcc Wigner crystal at low temperature, 
Carr^ did a calculation of the 3d Wigner crystal and de- 
termined the harmonic energy and the two electron ex- 
change integral. He did this assuming the ground state 
wave function is a product of single Gaussian orbitals 
and estimated the exchange frequency by calculating the 
exchange matrix element obtaining: 

Jnn ~ (1.6r;-^5 _ g_5^-i) exp(-1.55r]/2). (18) 

As shown in Fig. the slope of the exchange coefficient 
is roughly correct, while the prefactor is too small by a 
factor of about three. By including all effects of electron 
correlation, the tunnelling frequency is enhanced over 
what one gets with an uncorrelated wavefunction. More 
seriously, because of the two terms of opposite sign, he 
found an antiferromagnetic to ferromagnetic transition 
at Ts = 270. Carr's work is previous to the tunnelling 
theory of Thoulessi^ and Herringii, which asserts that 
purely pair exchanges must be antiferromagnetic, so a 
sign change like this can only come about from cancella- 
tion of exchanges of even and odd number of electrons. In 
Carr's calculation, there is no consideration of the possi- 
bility of more than two electron exchanging. The sign 
change comes from neglecting correlation between the 
electrons that are exchangingiii Given the strong cor- 
relations present at > 100, it is remarkable that the 
order of magnitude is reasonable. A reasonable estimate 
would not be obtained this way for the solid ^He because 
the hard core interaction makes the interaction matrix 
elements infinite, or nearly so. Carr also calculated the 
exchange for second neighbors; at Ts = 120 he obtained 
Jnnn/Jnn = 0.037, whilc wc obtained a ratio of 0.067. 



Now turning to comparison with other quantum sys- 
tems, similar behavior— in the spin Hamiltonian was 
found near melting between the second layer of ^He ab- 
sorbed on graphite and the 2dWC, but the similar be- 
havior does not seem to occur in 3d. It is suggested that 
the similarity is due to a common vacancy interstitial 
fluctuation mechanism, something that is also related to 
the melting of the quantum crystal^!, which is possibly 
second order or nearly so. 

This similarity is not found in 3d. The magnetic prop- 
erties of the 3dWC are quite different than bcc solid ^He. 
This comes about because pair exchanges are relatively 
more important in the Wigner crystal: we find in the 
3dWC that J*/ J„„ = 0.14 while in ^He the ratio is 0.41. 
This crucial enhancement of the pair exchange comes 
about because of the difference in the interaction strength 
at short distances: the helium-helium interaction has a 
much harder core. During a pair exchange, the particles 
have to pass by each other. As can be seen in Fig. |5]elec- 
trons can approach closer than helium atoms, and this 
allows them to pass by each other much more frequently. 
Exchange is a tunnelling process and the rate is very sen- 
sitive to the phase space in the transverse direction that 
the electrons have when they undergo exchange. 

The frustration between 2, 3 and 4-body exchanges 
leads, in solid helium, to a more complicated Neel or- 
dered ground state, the uudd phaseii. However, in the 
3d Wigner crystal, pair exchange dominates leading to a 
very stable antiferromagnetic phase. We find that the rel- 
ative frequency of 3- and 4-body exchanges is about the 
same in the two systems: Jt/Jp — 1.42 in ^He and 1.27 
in the 3dWC. In 3- and 4-body exchanges, the particles 
do not approach each other so closely, and the similarity 
in g{r) leads to similar ratios of exchange frequencies. 

Recently, it was found in DMC calculations'^ that the 
3d electron fiuid has a partially spin polarized ground 
state at low density. In fact, it undergoes a second order 
transition to a partially polarized state at = 50 ± 2. 
The polarization increases until it is fully polarized at 
freezing = 106. (Note, however, that because of the 
fermion sign problem, other types of ordering, such as su- 
perfiuidity, cannot be ruled out. In this respect, the liq- 
uid is much more difficult to treat than the crystal, since 
in the crystal, fermion effects are very small.) The region 
of stability of the spin polarized fermi liquid is shown in 
Fig. n It is curious that both the ffuid and crystal have 
magnetic ordering near the melting line. However, in the 
crystal, the magnetic ordering occurs at a temperature 
more than one thousand times lower than in the ffuid. 
As in the 2dWC, in the 3dWC the magnetic ordering 
temperature drops off very fast as density decreases. 

We hope to have provided a definitive result concern- 
ing the magnetic ordering of the 3dWC, a system which 
has provoked much speculation over the years. Because 
of the small energy scales, it will be an experimental chal- 
lenge to observe the magnetic ordering in the 3d Wigner 
crystal. 
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